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Abstract 

The properties of the wave function for the Modified Poschl- Teller (MPT) 
potential are outlined. The ladder operators are constructed directly from the 
wave functions without introducing any auxiliary parameter. It is shown that 
these operators are associated to the su(2) algebra. Analytical expressions for 
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the functions sinh(ax) and cosh ( ax ^ d_ are eva l ua ted from ladder operators. The 
expansions of the coordinate x and momentum p in terms of the su(2) generators 
^ ■ are presented. This analysis allows to establish an exact quantum- mechanical 

connection between the su(2) vibron model and the traditional descriptions of 
molecular vibrations. 
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1. Introduction 



During the past several decades, algebraic method has been applied to a wide variety 
of fields in both physics and chemistry. Systems displaying a dynamical symmetry can 
be solved with algebraic techniques [1-4]. Particularly, the Morse [5] and Poschl- Teller 
(PT) potentials [6, 7] represent two of the most studied anharmonic systems where these 
techniques have been used. Both of them are closely related with the 5*0(2, 1) [8-15] 
and SU(2) groups [3, 4, 16-18]. The latter has been used to describe the vibrational 
excitations of molecular systems, while the former is associated to the potential group 
approach. The relation between the SU(2) group and the Morse and PT systems 
can be directly established by means of a coordinate transformation applied to the 
radial equation of a 2D harmonic oscillator. The disadvantage of this approach need 
introduce an auxiliary variable to define the action of the group generators. It is 
possible, however, to obtain the ladder operators only in terms of the physical variable 
with the factorization method [16]. In this work, we shall present an explicit connection 
of the Modified Poschl- Teller (MPT) potential with the su(2) algebra. 

In the field of molecular physics, the PT potential was applied in the framework of 
the su(2) vibron model, where it was associated to the vibrational excitations of the 
molecular bending modes [19-21]. The basic idea behind the su(2) vibron model con- 
sists in expanding the Hamiltonian in terms of different sets of su(2) generators, whose 
number corresponds to the internal vibrational degrees of freedom of the molecule. 
Generally, the stretching degrees of freedom are associated to a Morse potential but 
the bending coordinates can be identified with either the Morse or PT potentials. This 
model was empirically applied without establishing its connection with the traditional 
approach in the configuration space. Consequently, in such a scheme the relation be- 
tween the spectroscopic parameters and the molecular structure constants was lacking. 
Concerning the degrees of freedom associated to the Morse potential , this problem has 
been recently solved by establishing the relation between the matrix elements of two 
interacting Morse oscillators and the corresponding SU(2) matrix elements. The su(2) 
vibron model was shown to correspond to an approximation where only the dominant 



2 



Av = ±1 interaction between Morse oscillators is taken into account [17]. Moreover, 
a procedure to extend the algebraic approach in order to systematically approach the 
Morse oscillators results was proposed [18]. This procedure is based on the expansion of 
the coordinate x and momentum p in terms of the su(2) generators; an expansion that 
was obtained by comparison of the matrix elements in both spaces. In this work we 
shall obtain the corresponding expansion of x and p for the MPT potential by inverting 
the expressions for the ladder operators. This analysis provides us for an approach to 
approximate the matrix elements of interacting PT oscillators in the framework of the 
su(2) model. 

This paper is organized as follows. The properties of the MPT potential are outlined 
in Section 2. In Section 3, the associated ladder operators are constructed directly 
from the wave functions with the factorization method [22-24] and their connection 
with the su(2) algebra is established . In addition, the matrix elements for certain 
relevant functions of the coordinate x and momentum p are calculated from the ladder 
operators. In Section 4 the expansions of the coordinate x and momentum p in terms of 
the su(2) generators are presented , including an analysis of two coupled PT oscillators 
system. Some concluding remarks will be given in Section 5. 



2. The properties of the MPT potential 

We start by presenting the solutions of the MPT problem [25]. The MPT potential 
as described in Ref. [26] can be written as 

where D is the depth of the well and a is related with the range of the potential, while 
x gives the relative distance from the equilibrium position. The Schrodinger equation 
associated to this potential is given by 

d 2 m(x) 2a / D \ T „ . . 
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where ji is the reduced mass of the molecule and q is related with the depth of the 
potential. The solutions of Eq. (2) can be written as 

2fc2 

W n (u) = JV*(1 - u 2 )iC q n +h - n (u), E n = -%—{q - n) 2 , n = 0, 1, 2, ... (3) 



2// 



where 



where 



" /4jE ,g(g + l) = ^,g = i(-l + 2A;),A; = J] + = 2k = 2q + l, (4) 



2 ft 2 ' yw ; a^ 2 ' H 2 V v ' V4 a 2 h 2 ' 



u = tanh(a:r), e = q — n > 0, 



an!(g — n — \)\{2q — 2n)\ 

(5) 



\ ?r3(g-n-l)!(2g-n)! 
. Here z/ has been introduced because of its relevance for the identification of the 

ladder operators with the su(2) algebra, as will be shown in the next section. On 

the other hand, the number of bound states is determined by the dissociation limit 

e = q — n = 0. We should note that for q integer the state associated with null energy 

is not normalizable. In this case the last bounded state corresponds to q — n — 1. We 

thus have that n max — q — 1 = {v — 3)/2. 



3. The construction of the ladder operators 

In this section we address the problem of finding ladder operators with the fac- 
torization method , namely, we intend to find the differential operators V with the 
following property 

V±KW=P±f q n ±l(u). (6) 

Specifically, we look for operators with the following form 

V± = A ± {u)^ + B ± {u) ) (7) 

where we remark that these operators only depend on the physical variable u. 

The operators (14) can be found by obtaining the action of the differential operator 
on the MPT wave functions. Therefore the following formula (see Eq. 8. 935 in 

[28]) 



^ = 2A<#m (8) 
together with Eq. (3), allows to obtain 

and introducing the explicit form of the normalization constant, equation (9) becomes 



v / T^ (± + ^-A) f-^^niu) = y/nfr-n + lW-M, (10) 
from which we can define the annihilation operator P_ as 



^ (id 



du 1 — -u 2 / V 9 — n 



or in terms of z/ defined in (5) 



'-^(s + i^'Jf^ (12) 

where in order to simplify the notation we have taken into account that 2e = v— 2n— 1 = 
2g — 2n. The action of the operator (12) on the wave functions is then given by 

P-KM = p-K-iM, (is) 

where 



p_ = \jn{y — n). (14) 

As we can see, this operator annihilates the ground state ^(u), as expected from a 
lowering operator. 

We now proceed to find the corresponding creation operator P + . To this end, we 
consider the formula (see (10. 43) of [29]) 

2(\-l)(2\-l)xC*(x) = 4\(\-l)(l-x 2 )C*+l(x) + (2\ + n-l)(n+l)C*-l(x). (15) 
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This recurrence relation can be used together with Eq. (8) to obtain 
dn(u) _ u(g - n) (n + l)(2g - n) iVg 

By using the explicit form of the normalization constant (5), equation (16) becomes 



_ rii2 ( d u(q-n) 



K z J ] ^ u ) = + - nW n+1 (u). (17) 



du 1 — u 2 J \ (g — n) 
Likewise, in terms of the variable z/, we can thus define the creation operator P + as 



,2 



d u \ e — 1 



r/// 1 if 2 j v 

with the following effect on the wave functions 

P + ^(«)=p + *^ +1 («), (19) 

where 



p + = > /(n + l)(i/-n-l). (20) 

From Eqs. (8) and (17) we see the importance of normalization factors. Since P + is a 
raising operator it is expected to annihilate the last bounded state. Indeed, for such 
state e = 1 and the square root in (18) makes the operator vanish. 

We now establish the algebra associated with the operators P±. Based on the Eqs. 
(12) and (19), we calculate the commutator [P_,P + ]: 

[P + ,P_]*» = 2po^(«) (21) 
where we have introduced the eigenvalue 

Po = ~ ( - n) . (22) 
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We can thus define the operator 
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P =n- V -±. (23) 
The operators P± ; o satisfy the commutation relations 

[P+,P_]=2P , [P ,P-] = -P-, [Po,P + ] = P+, (24) 

which correspond to the su(2) algebra. This result is consistent with the description 
of a finite discrete spectrum, in accordance with previous algebraic descriptions of the 
bounded states of the Poschl- Teller potential [4, 8] . The number of bosons N is related 
with v by ./V = v — 1, as we deduce from the Casimir operator 



C WJu) 



K(u)=j(j + l)K(u), (25) 



where j, the label of the irreducible representations of su(2), is given by 

From the commutation relations (24), we know that P is the projection of the angular 
momentum m, and consequently 

v- 1 

n - — = m. (27) 

The ground state thus corresponds to m — —j, while the maximum number of quanta 
n max = and consequently m max \ nmax = —1, in accordance with the constraint 
condition e = q — n = 1 for the last bounded state. The MPT wave functions are 
thus associated to one branch (in this case torn < —1) of the su(2) representations, as 
expected. Finally we should notice that in terms of the su{2) algebra, the Hamiltonian 
acquires the simple form 

H = -^Pl (28) 



where 



h(3 2 u 

w — — • 
2/i 
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while for the wave functions 



(29) 



where the normalization constant is obtained through the commutation relations (24), 
and turns out to be 



AC = 



(u-n- 1)! 



(30) 



\| n\(u-l)\ ' 

For other calculations one can obtain the following expressions in terms of the 
raising and lowering operators P± 



u 



= - P- 



e(e + l) 



+ ^4 



c(c-l). 



(31) 



(32) 



du 2\~ V e+ 1 ~ T Y e ~ 
where it has to be understood that for the last bounded state (e = 1) the raising 
operator vanishes. On the other hand, we remark that the variable e is to be considered 
as an n dependent operator. This is the price to pay for having ladder operators in 
terms of the physical variable u, which contrasts with the 2D description, where no 
number operators are involved but an auxiliary variable must be included [13]. Using 
Eqs. (13) and (19) and considering the constraint condition 2e = v — 2n — 1, we can 
thus calculate the matrix elements of these functions as 



<*n'Wl7=?l*M> = <K'(x)\s™Hax)\W n (x)> 



+ 



n{y — n) 



\ {y-2n- l)(i/-2n + l) 



"n'.ji— 1 



(n + l)(v -n - 1) 
\\ (i/-2n-l)(i/-2n-3) 



n',n+l 



(33) 
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<K^)\VT^£\K(n)> = 



I cosh(a:r) d 
a dx 



K > M 



n{y — n)(u — 2n — 1) 
2^ {y-2n + 1) 



(34) 



{n + l)(v-n - \)(y-2n- 1) 
(i/ - 2n - 3) 



+ 1 



5. Connection with the su(2) vibron model 

In the previous sections it is showed that the su(2) algebra can be identified as 
the dynamical group for the MPT potential. Introducing the following renormalized 
expressions 

St = ^, b = ^==, k = —. (35) 
we can obtain the simplified MPT wave functions 



CM 



' l! (^"W. (36) 



\| n!(i/-l)! 

which implies that he wave functions can be expressed in terms of a repeatedly action of 
the creation operator w on the ground state. On the other hand, it was well known the 
connection of the PT potential with the SU (2) group through a 2D harmonic oscillator. 
This connection allows to associate the PT potential to certain bending degrees of 
freedom (e.g. the out of plane local modes) in the framework of the su(2) vibron 
model. This model consists in expanding the Hamiltonian in terms of su(2) generators 
associated to local oscillators (internal coordinates). It is not obvious, however, the 
relation of such generators with the creation and annihilation operators obtained in 
Section 3. The aim of this section is to establish this connection. 

Let us consider two interacting oscillators in the framework of the su(2) vibron 
model. According to this approach the Hamiltonian has the form [18] 

H su{ 2) = \^ + Kb\) + \hu (b\b 2 + M|) . (37) 

i=i 



The first term corresponds to the independent oscillators identified with either the 
Morse or PT potentials. In this study, we shall consider PT potential. The parameter 
uo and the boson number N are related to the PT spectroscopic parameters through 



[18] 



where 



uj e = Tiu) = Huq 



N 



hu 



X e UJ e — Ar 

N 



(38) 



_ u e a 2 h 1 2uD 



2*2 



a 2 h 



x e u e 



The diagonal terms generate PT (Morse)-like spectra, while the interaction term 

iC (2) = Afta* (b\b 2 + hbl) , 
gives rise to the matrix elements 

< N,m+1, N,n 2 -l\HZ {2) \Nn u Nn 2 >= XhuJoJr^jH + 1) 



(39) 



(40) 



(n 2 - 1) 



N 



1-^ 

N. 



(41) 

where n« stands for the number of quanta of the i-th oscillator. In the su(2) vibron 
model the local functions are denoted by \N,n >. We note that corrections of order 
1/N appear. When N — > oo the matrix elements (41) reduce to 

lim < m + l,n 2 - l\H-nt^\ni,n 2 >= XHloqJ n 2 {ni + 1) , (42) 

as expected, since in this limit all the harmonic results are recovered [18]. In order to 
interpret the correction appearing in (42), we shall analyze two PT oscillators from the 
point of view of the dynamical group identified in Section 4. 

The Hamiltonian of two identical MPT oscillators in configuration space is given 

by 



1 2 2 2 1 



+ A \ — — + jjiw r x r 2 I . (43) 



The first two terms correspond to two PT independent oscillators, while the third is a 
particular interaction chosen to correspond in the limit N — > oo to the interaction 



Htt = ^ (a\a 2 + MJ) 



(44) 
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for two harmonic oscillators. In configuration space Eq.(44) takes the equivalent form 

H£ = \^+ijuJ 2 x 1 x^ , (45) 

with uj given by 

^ = \lfrr9rv (46) 

Note that the matrix elements of the interaction (44) coincide with the matrix elements 
in (42). 

We now address the problem of computing the non diagonal contributions to the 
interaction term 

H int = A I — — + fjw x x x 2 I . (47) 
To this end we introduce the bosonic operators 




In terms of these operators the interaction (47) acquires the simple form 

H^ nt PT = hu\(Ac2 + cA) ■ (49) 

In order to establish the relation between this term and the interaction associated to 
the su(2) vibron model given by (40), we proceed to express (49) in terms of the su(2) 
ladder operators obtained in Section 3. This task can be achieved by extracting the 
coordinate x and momentum p from (31) and (32), which implicitly are given by 



sinh ax = \ (V^TT) + *V^i)) ' (50) 



coshfarr) d 1 / f / ve ? + / ve \ . . 

b \hrr- b \ — -r • ( 51 ) 



a dx 2 y y e + 1 V e — 1 
If we now take into account the Taylor series 



arcsinh(x) — x — — + + (52) 
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sech(x) = 1 - — + — + 



(53) 



keeping up to the second terms, we obtain from (50) an (51) the following expansions 

ax = ]^{b ] f n + bg n ) - ^(b ] f n b ] f n b ] f n + o f f n b ] f n bg n + b ] f n bg n b ] f n + b ] f„bg n bg n 
+ bg n tff n tff n + bg n Pf n bg n + bg n bg n P f n + bg n bg n bg n ) + ... 



(54) 



p = l l^(jjh n - bq n ) + l -^(b ] f n b ] j r pK - b ] fjff n bq n + b ] f n bg n tih n - o f f n bg n bq r , 
2 lb 

+ bg n Pf n b ] h n - bg n Pf n bq n + bg n bg n b ] h n + bg n bg n bq n ) + 



(55) 



where we have introduced the notations 



fn 



V 



6(6-1)' 



V 



9n 



e(e + l) 



(56) 



h„, = 



V e-1' w "V e + 1 ' ^ 
as well as p = —ih-^. We recall that the n dependence if these variables are given 
through 2e = v — 2n — 1. Note that because of the symmetry of the MPT potential 
the expansions of x and p involve only odd powers of the operators b^(b); x and p are 
odd functions under inversion. If we introduced the expressions (54) and (55) in the 
definitions (48), we would obtain an expansion of c*(c) in terms of the P(b) operators. 
As a first approximation, however, we shall consider contributions up to linear terms 
in (54) and (55). If this is the case the operators (48) take the approximate form 

ha 




2a 



(tff n + bg n ) 



+ 



■{b ] h n - bq n ) 



(58a) 



c ~ 




2a 



{Pfn + bg n ) 



2Tiujji 



ha 



(tfh n - bq n ) 



(586) 
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which can be rearranged in the form 



C f ~ P Z n + 6 Cn, 



(59a) 



c ~ 6 f Cn + b z n , 
where we have introduced the functions 
1 



1 



1 - (2n + l)/u 

2 I V (1 - (2n + - (2n + 3)/z/) ' \\ 1 - (2n + 3)/v J ' 



+ 



1 - (2n + \)/v 



2 \ V (1 - (2n + - (2n - \ 1 - (2ra - 1)/ 

It is clear that these functions have the harmonic limit 

lim z n = 1, lim Cn = 0, 



(596) 



(60) 



(61) 



(62) 



although for small number of quanta this limit represents a good approximation. Only 
when n takes its highest value this approximation fails. Hence we have the following 
limit 

lim $ = b\ lim c — b. (63) 

Let us now return to the two oscillators system. Considering the approximation (59), 
we obtain 



\c 2 + cicj) ^ (6162 + SiSJ)i(m, n 2 ) + (6I&1 + SiS2)S(m, n 2 ) . 



(64) 



where 



A(m,n 2 ) = 4 X ^n 2 + Cm Cn 2 , (65) 
5(ni,n 2 ) = z m 

C«2 ^"2 Cm- (66) 

From the discussion for one oscillator, a reasonable approximation for small ri\ and n 2 
is the following 

A(m,n 2 )~l, 5(m,n 2 )~0, (67) 

which implies 

#^ PT = hw\(c\c 2 + ci^) - «wA(6 t 1 6 2 + 61SJ) . (68) 
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This expression corresponds essentially to the su(2) vibron interaction (40). We thus 
conclude that the su(2) vibron model very nearly corresponds to taking the dominant 
An = ±1 interaction between coupled MPT oscillators, as long as we exclude the 
neighborhood of the dissociation energy. 

This analysis suggests the following two steps to improve the the su(2) vibron 
model. First we can consider the operators (59) instead of the single generators W and 
b. The next step would include the next terms in the expansions (54) and (55). Both 
steps imply the breaking of the polyad P = ri\ + n 2 in a specific way, specially for high 
number of quanta. This approach may thus be useful for the descriptions of vibrational 
excitations near the dissociation limit. 



5. Concluding remarks 

In this work we have outlined the properties of the MPT wave functions and estab- 
lished the ladder operators with the factorization method. The realization obtained 
for the creation and annihilation operator, given in terms of the physical variable u, 
has been identified with the su(2) algebra. This result is in accordance with previous 
analysis in which the bounded region of the MPT potential is associated with an su(2) 
algebra. The identification of the ladder operators as generators of the SU(2) group 
allows to express the MPT wave functions in a simple closed form in terms of the 
action of ¥ on the ground state. The matrix elements for the functions sinh(a::r) and 
cosh(ax) d_ j^yg been analytically obtained from the ladder operators P±$. This method 
represents a simple and elegant approach to obtain these matrix elements. 

We have also analyzed the theoretical relationship between the su(2) vibron model 
and the coupled MPT system. It has been shown that the standard su(2) approach 
to molecular vibrational excitations (vibron model) is equivalent to considering MPT 
oscillators via the dominant An = ±1 selection rule. This approximation arises nat- 
urally from an expansion of the PT operators $(c) in terms of the su(2) generators 
W(b), which in turn allows to establish an extended su(2) model when the next terms 
in the expansion are taken into account. The extended su(2) model may provide an 
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approach to the vibrational description of the high energy region of the spectra, where 
polyad is breaking and localization takes over the simple normal behavior at the low 
end of the spectrum. 
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